Cold guided beams of water isotopologs 
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Electrostatic velocity filtering and guiding is an established technique to produce high fluxes of 
cold polar molecules. In this paper we clarify different aspects of this technique by comparing ex- 
periments to detailed calculations. In the experiment, we produce cold guided beams of the three 
water isotopologs H2O, D2O and HDO. Their different rotational constants and orientations of 
electric dipole moments lead to remarkably different Stark shift properties, despite the molecules 
being very similar in a chemical sense. Therefore, the signals of the guided water isotopologs differ 
on an absolute scale and also exhibit characteristic electrode voltage dependencies. We flnd excel- 
lent agreement between the relative guided fractions and voltage dependencies of the investigated 
isotopologs and predictions made by our theoretical model of electrostatic velocity filtering. 

PACS numbers: 37.10.Mn, 37.10.Pq, 32.60.-l-i 
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INTRODUCTION AND MOTIVATION 



Cold polar molecules offer fascinating perspectives for 
research, e.g. in cold and ultracold chemistry, preci- 
sion measurements for tests of fundamental symmetries, 
and quantum information (see, e.g., the special issue 
on cold polar molecules Since molecules are in 

general inaccessible to direct laser cooling, new meth- 
ods are needed. Indirect methods such as association 
of molecules from ultracold atomic ensembles by pho- 
toassociation [1, d, [1] or using magnetic Fcshbach res- 
onances 0, [1,13| have the advantage of producing po- 
lar molecules directly at ultralow translational temper- 
atures. Although the molecules are produced in highly 
excited vibrational states, it has been shown that it is 
possible to transfer them down the vibrational ladder to 
create stable ultracold molecules 0, H, [13, UH. How- 
ever, these techniques are limited to a few species which 
can be forged together from laser-coolable atoms. Direct 
methods, which can be applied to naturally occurring 
polar molecules, include buffer-gas cooling [l^, electric 
111, [13, [Hi and optical [TgI . [TtI Stark deceleration, mag- 
netic deceleration [Tsl . Il9|. c ollisions with counterprop- 
agating moving surfaces [23| or collision partners [2]| . 
rotating nozzles [l^ , velocity filtering of large molecules 
by rotating mechanical filters and velocity filtering 
by an electrostatic quadrupole guide [U, [11] . 

Electrostatic velocity filtering offers the advantage of 
being a simple technique delivering continuous beams of 
high flux; guided beams of more than ipio molecules/s 
have been produced for ND3 and H2CO (gj, i25|]. This 
electrostatic guiding and filtering technique is well suit- 
able for collision experiments [26| . since it produces a 
high continuous flux with only few rotational states con- 
tributing [23, [ll]. ft has been combined in a natural 
way with buff er-g as cooling to increase the purity of the 
guided beam [29|. Furthermore, it is applicable to many 
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molecular species, as long as they have populated low- 
field-seeking (Ifs) states in the thermal source. For high- 
field-sceking (hfs) states, velocity filtering in the guide is 
also possible by changing from a static quadrupole field 
to a time-dependent dipole field (30j . 

The filtering by an electric guide strongly depends on 
the Stark shift properties of the molecules used. In gen- 
eral, dissimilar molecules not only have different Stark 
shifts, but also different physical and chemical proper- 
ties, which play an important role in the source and in 
the detection process. Therefore, it is difficult in gen- 
eral to compare such measurements. In this paper we 
present experiments performed with the three water iso- 
topologs H2O, D2O and HDO. These molecules have 
comparable masses and chemical properties but differ in 
their rotational constants and dipole moments. This al- 
lows the study of the velocity filtering process without 
large uncontrollable systematic effects. Although these 
isotopologs seem very similar, they show surprisingly dif- 
ferent behavior when exposed to external electric fields. 
Their different Stark shift properties are experimentally 
revealed as characteristic dependencies of their detector 
signals and velocity distributions on the applied electrode 
voltages. The experimental work is accompanied by a 
theoretical description of the filtering process. We find a 
good agreement between the predictions of our calcula- 
tions and the experiments, concluding that the filtering 
process is well described by the presented model. 

The paper is organized as follows: In Section |TT] we 
review the filtering process in an electric guide for polar 
molecules. In Section Hill we present the theory for Stark 
shifts of asymmetric rotor molecules, which is applied to 
the water isotopologs H2O, D2O and HDO. In Sect ion HVl 
it is shown how the calculated Stark shifts of H2O, D2O 
and HDO and the described filtering properties of the 
guide can be used to make predictions of the guided flux. 
The experimental setup used for the guiding experiments 
is described in Section |Vl From detailed measurements 
presented in Section IVII we determine relative detector 
signals and velocity distributions of cold and slow water 
molecules. 
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II. VELOCITY FILTERING 

To calculate the relative guided fluxes of the differ- 
ent water isotopologs, it is necessary to review some ba- 
sics of velocity filtering by an electric guide. For guid- 
ing of low-ficld-sccking molecules, the guide electrodes 
are charged to positive and negative high voltages in a 
quadrupolar configuration (as shown in Fig. [8] in Section 
|V| . This creates an electric field minimum in the center, 
surrounded by a linearly increasing electric field, up to 
a certain maximum trapping field. In this electric field 
minimum molecules in Ifs states can be trapped in trans- 
verse direction and guided. 

Molecules from a thermal reservoir at temperature T 
are injected into the guide. Their velocities are described 
by a three-dimensional Maxwell-Boltzmann velocity dis- 
tribution 



-u^ exp(— 



f{v)dv = - 

and by a one-dimensional velocity distribution 
1 



(la) 



■cxp(-u^ /a )dv^,y^^ (lb) 



with most probable velocity a ~ yj2kBT jm, veloc- 
ity components Vi {i = x,y,z) and total velocity v = 

+ Vy + v1 . The velocity distribution of molecules 

coming from the nozzle is given by 



P{vz)dvz = -^Vz exp{-vl/a'^)dvz 



(2) 



where the z-direction is defined to be oriented along the 
guide. The factor enters since one is now considering 
the velocity distribution in the flux out of the nozzle and 
not in a fixed volume. 

The guided fraction of molecules can be calculated as 
the part of molecules injected into the guide with trans- 
verse and longitudinal velocities below certain trans- 
verse and longitudinal cut-off velocities, which depend 
on the Stark shift of the molecules and the properties 
of the guide. For a molecule with a given Stark shift 
AW'^{Emax) at the maximum trapping field Emax, a 
maximum transverse velocity Vmax = ^/2AW^/m ex- 
ists. If the transverse velocity v± of the molecule ex- 
ceeds Vmax, it is lost from the guide. In the follow- 
ing it is assumed that the molecule is injected in the 
center of the guide. If molecules enter the guide off- 
center they have acquired already some potential en- 
ergy which reduces the maximum trappable transverse 
velocity. Furthermore we assume mixing of the trans- 
verse degrees of freedom, whereas longitudinal and trans- 
verse degrees of freedom are not coupled. In a full nu- 
merical analysis this assumption is dropped and yields 
only small modifications to the guided fiux. Filtering 
on longitudinal velocity vi is realized by bending the 
guide. In the bend, molecules with a velocity below the 



maximum longitudinal guidable velocity vi^max remain 
trapped. The longitudinal and radial maximal veloci- 
ties are connected by the centrifugal force acting on the 
molecules in the bend, leading to vi^max oc Vmax- More 
specific, the guide has a free inner radius r, at which 
the maximum trapping field is reached. Since the trap- 
ping potential increases linearly in radial direction, the 
maximum restoring force acting on a molecule is given 
by Fr^max = 'AW ^ (Emax) / f ■ The maximum longitudinal 
velocity is then obtained by equating the centrifugal force 
in the bend of radius R and the restoring force, resulting 
in vi^rnax = \J AT^ {Emax)R/r m = y/R/2rv,nax- Note 
that for increasing longitudinal velocity the maximum 
trappable transverse velocity decreases. As basically ev- 
ery particle is guided if vi < vi^nax and v± < Vmax, this 
results in higher efficiencies as compared to filtering by, 
e.g., rotating filter wheels and apertures. 

To calculate the guided flux, the velocity distributions 
of molecules injected into the guide are integrated to the 
maximum trappable velocity. When performing the in- 
tegration in the limit of small cut-off velocities Vmax and 
Vi^max compared to the thermal velocity a, the expo- 
nential cxp(— u^/a^) in the thermal velocity distributions 
can be replaced by 1. The guided flux $ of a molecular 
state with a Stark energy AW* is then given by 

$ = 



.f{Vx)fiVy)PiV;:)dVxdVydVz 

Vx— — f^max Vy— — Vmax ^5—0 

««L,(x(AW)2, (3) 

where we have neglected the dependence of the trans- 
verse cut-off velocity on the longitudinal velocity, and <& 
is normalized to the flux out of the nozzle. The guided 
flux can hence be described by a function / which gives 
the fraction of guidable molecules for a given electric field, 
$ = /(AVF") (X (AW)2. The dependence of the guid- 
able fraction given by /{AW) and therefore the depen- 
dence of the guided flux on the molecular Stark shift al- 
lows us to infer characteristic Stark shift properties from 
measurements of the guided flux. When varying the elec- 
trode voltage and hence the guiding electric flcld, the flux 
will change depending on the molecules' Stark shift. As 
will be shown in Section [TTTl molecules can exhibit linear 
or quadratic Stark shifts, depending on the molecular 
properties. In Section IIVI the influence of these Stark 
shift behavior on the guided flux will be discussed in de- 
tail. In short, the flux scales quadratically with electrode 
voltage for molecules with a linear Stark effect and quar- 
tic for molecules with a quadratic Stark effect, as can be 
seen from Eq. Q and as already described in [2^, [2^ . 



III. STARK SHIFTS 

As shown in Section [Til the guided flux of molecules 
in a certain rotational state is determined by the Stark 
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D2O 


Rotational 
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27.79 


23.48 


15.39 


constants 
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14.50 


9.13 


7.26 


[cm-i] 
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9.96 


6.40 


4.85 


Dipole moment 







0.66 





components 


fJ-B 


1.94 


1.73 


1.87 


[Db] 


fJ-c 












TABLE I: Rotational constants and components of the elec- 
tric dipole moment along the principal axes for the different 
water isotopologs H2O, D2O and HDO [3^. |33|. 




FIG. 1; (Color online). Orientation of the rotational axes in 
D2O. The electric dipole moment is oriented along the 6-axis. 

shift of that state. The total flux is then given by the 
sum over the contributions of all states thermally popu- 
lated in the source. To make a theoretical prediction for 
guided fluxes of the different water isotopologs H2O, D2O 
and HDO, their Stark shifts are calculated by numerical 
diagonalization of the asymmetric rotor Hamiltonian in 
the presence of an external electric field, following the 
procedure given by Hain et al. [3l|. For completeness, 
the main steps of the calculation are summarized in Sec- 
tion IIII Al A detailed understanding of this theoretical 
part of the paper is not crucial for the discussions follow- 
ing. Section IIII Bl shows how the molecular properties 
manifest themselves in the Stark shifts, and illustrates 
this with the example of the different water isotopologs. 

A. Calculation of Stark Shifts 

The molecular properties relevant for the calculations 
are listed in Table [J A, B and C are the rotational con- 
stants along the principal axes, where A (C) is oriented 
along the axis with largest (smallest) rotational constant, 
A > B > C . The orientation of the rotational axes in 
D2O is illustrated in Fig.[TJ 

The molecular Hamiltonian in the absence of external 
electric fields can be written as 

Hrot^\{A + C)3-' + \^{A-C)H{k), (4a) 



with the reduced Hamiltonian 

H{k) =31 + Kil - 31 (4b) 

The constant k = {2B -A- C)/{A - C) is the so called 
asymmetry parameter, taking on the value —1 in the limit 
of the prolate symmetric top and -f 1 in the limit of the 
oblate symmetric top. For calculations this Hamiltonian 
is expressed in the symmetric rotor basis. The matrix 
elements of the reduced Hamiltonian H{k) [Eq. (j4bp ] in 
the symmetric top basis {| JA'M)} are given by 

{JKM\ H{k) \JKM) = F [J{J + 1) - + GK^, 

(5a) 

{JJ<±2,M\H{n)\JKM) =i^[/(J,i^±l)]^/^ (5b) 
with 

f{J,K ±1) = -^[J{J + 1) - K{K ±1)] 

y.[J{J +l)-{K±l){K±2)] (5c) 

and F, G, H supplied in Table Ivl and I VII in the appendix. 
Only symmetric rotor states with /S.K = 0, ±2 are cou- 
pled by the asymmetric rotor Hamiltonian as is seen from 
Eqs. dni). The eigenstates Ajru and energies WjrM of 
the asymmetric rotor can be found by diagonalization of 
the full asymmetric rotor Hamiltonian Eq. (|4ap . yielding 

WjrM = \{A + C) J(J -M) + 1(A - C)WjrM{K) (6) 

and 

AjrM^Y.''K^''^JKM. (7) 
K 

The eigenstates Ajtm are expressed as linear superpo- 
sition of symmetric rotor wave functions ^ jkm- Note 
that the total angular momentum quantum number J 
and its projection on a space-fixed axis M are still good 
quantum numbers in the field-free asymmetric rotor, in 
contrast to K . This can already be seen from the asym- 
metric rotor Hamiltonian which does not depend on M 
and which conimutates with J^. The pseudo quantum 
number r = Ka — Kc is used to label the asymmetric 
rotor states in ascending order in energy, r is directly re- 
lated to the quantum numbers Ka and K,. in the prolate 
and oblate limiting case of the symmetric top. From this 
the symmetry properties of the asymmetric rotor states, 
represented by the symmetry species of the four group 
D2 nil (sometimes also referred to as V{a,h,c)), can be 
derived [H,!!^. 

An external field lifts the degeneracy between the M- 
sublevcls of a state | J, r) . Different J states are also cou- 
pled now, leaving M as the only good quantum number. 
In general, an asymmetric top molecule can possess com- 
ponents of its dipole moment /i along all three principal 
axes in the body-fixed frame , pi = l^g^g^ {q — c). 
As can be seen from Table HI H2O and D2O have a dipole 
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moment component only along the 6-axis (see Fig. [T]), 
whereas in HDO there are components along the a- and 
6-axis. This is important for their Stark shift properties, 
since different components promote couplings between 
rotational energy levels of different symmetry species. 
We define the electric field to be directed along the Z- 
axis in the space-fixed frame {F = X, Y, Z), E = Ez&z- 
Then, the interaction Hamiltonian is given by 



-M=0 



= E, 



(8) 



where the direction cosines ^Pg connect the space- fixed 
to the molecule-fixed frame. Since the direction cosines 
are tabulated for symmetric rotor wave functions only 
(see Table IVIII in Appendix „ note there are some mis- 
prints in the <^pg tabulated in [3l| which are corrected 
here.), their matrix elements with respect to the asym- 
metric rotor states have to be constructed. This can be 
done using the expansion of the asymmetric rotor wave 
functions in terms of symmetric rotor states and results 
in 

{jTM\'^Zg\J'T'M') = 

{J\'^Zg\J') (JM|$Z<,|J'A/') 5mM' 



M^J't'M' 



a 



K' 



\JK\^Zg\J'K'). (9) 



KK' 



Note that Hs only couples states with AAf = 0, due to 
the choice of E along the Z-dircction, E = Ez^z- I3y 
diagonalization of the total Hamiltonian Hrot + Hs{Ez) 
the eigenstates and eigenenergies of the asymmetric ro- 
tor in the presence of an external field can be calculated. 
The Stark shift AW'^{Ez) is then simply given by the dif- 
ference between the total energy in the external electric 
field and the zero field energy. 

B. Discussion of Stark Shifts 

The energies of the rotational states of the different 
water isotopologs in an external electric field, calculated 
by the procedure described in Section Iffl Al arc shown 
in Fig. [5] and [31 Several features are directly evident 
when comparing the different isotopologs. First of all, 
the calculations predict quadratic Stark shifts for H2O 
and D2O, while for HDO linear Stark shifts are found as 
well. This is caused by the different orientations of the 
electric dipole moments in the molecules with respect to 
the main axes of the molecule. In H2O and D2O the 
dipole moment is only oriented along the 6-axis. In HDO 
there is additionally a component along the a-axis. This 
is important since in a more prolate asymmetric rotor, 
as is the case for the different water isotopologs (H2O 
K = -0.49, D2O K = -0.54, HDO n = -0.68), the states 
with the same Ka quantum number are near degenerate 
for increasing J and Ka quantum numbers. It is exactly 
these states which are coupled by a dipole moment along 
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FIG. 2: (Color online). Energy of rotational states in an 
external electric field for the water isotopologs H2O, D2O and 
HDO. The same vertical scale is used throughout the figure. 
States with large contributions in the guided beam are plotted 
with thick red lines and marked by a star. 



the a-axis in an external electric field. Hence, once the 
coupling between these states due to the electric field 
becomes larger than their asymmetric rotor splitting the 
Stark shift becomes linear. This behavior is illustrated 
in Fig. [3] for the lowest energy rotational states of HDO 
exhibiting a mainly linear Stark shift. For H2O and D2O 
the situation is completely different, here the Stark shifts 
stay quadratic and even become smaller with increasing 
rotational energy as can be seen from Fig. [D 
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\J,Ka,K,) 




0.25 



Electric field (kV/cm) 

FIG. 3: (Color online). Lowest energy rotational states with 
linear Stark shifts of HDO. The same vertical scale is used 
throughout the figure. States with large contributions in the 
guided beam are plotted with thick red lines and marked by a 
star. In this figure quantum numbers Ka and Kc are used to 
indicate the close connection to the prolate symmetric rotor 
where these states are degenerate, giving rise to linear Stark 
shifts. 



A second observation which can be made from Figs. Hi- 
ll] is that the magnitude of the Stark shifts of the iso- 
topologs also differ. HDO exhibits the largest Stark 
shifts, as is to be expected because of their linear char- 
acter. Besides, the Stark shifts of H2O are significantly 
smaller than those of D2O, although both molecules have 
similar size dipolc moments and hence similar couplings 
between rotational states. The reason for this is also 
found in the rotational constants. Since D2O has smaller 
rotational constants than H2O, the energy levels are 
closer together, leading to larger Stark shifts for compa- 
rable couplings, as can be seen from perturbation theory. 

The third observable feature is that the Stark shifts 
also show a different behavior with increasing rota- 
tional energy for the different isotopologs. In HDO, 
the states | J, t, Af) = | J, J, J) exhibit the maximum Stark 
shifts which approach a constant value given by the max- 
imum possible projection of the dipole moment compo- 
nent along the a-axis on the electric field axis as in- 
dicated in Fig. H The Stark shifts in H2O and D2O de- 
crease with increasing rotational energy, as was alr eady 
observed and described for D2O in Rieger et al. [28j . 
With increasing J and K quantum numbers the spac- 




Rotational energy (cm ') 

FIG. 4: (Color online). Stark shifts of the different wa- 
ter isotopologs H2O, D2O and HDO at an electric field of 
130kV/cm. The solid line indicates the Boltzmann factor 
for a source temperature of 293 K. Only Ifs states are shown. 
Note the different vertical scales. The different M states of a 
rotational state \J,t,M) all have the same rotational energy 
but different Stark shifts, thus forming a vertical sequence. 
Note that in the case of HDO the Stark shifts approach a 
saturation value indicated by the dashed line for increasing 
rotational energy, corresponding to the full alignment of the 
dipole moment component along the a-axis on the electric 
field axis. In the case of the fo-type rotors H2O and D2O the 
Stark shifts decrease with increasing rotational energy. 



ing between the rotational energy levels increases, which 
reduces the Stark shifts. 

More generally, the Stark shift properties discussed 
above can be understood for an arbitrary asymmetric top 
molecule by examining its rotational energy level struc- 
ture, the symmetry properties of its energy levels, and 
the couplings brought about by the different components 
of the molecular electric dipole moment. These couplings 
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Symmetry Species: A B, 

a) c) 




FIG. 5: (Color online). Couplings due to the different dipole moment components in the asymmetric rotor. Shown are the 
energy levels of the asymmetric rotor as a function of the asymmetry parameter k, for different J quantum numbers. The 
symmetry of the rotational states, represented by the symmetry species of the four group D2, is indicated by the dashing of 
the curves. A dipole moment fia along the a-axis mainly couples states which are degenerate in the prolate symmetric top, i.e. 
K = —1, while a dipole moment along the c-axis mainly couples states degenerate in the oblate symmetric top, i.e. k = 1. 
This near degeneracy gives rise to linear Stark shifts for these states. In contrast, states which are coupled by a dipole moment 
/ii, along the 6-axis, are always separated by an intermediate energy level, therefore never giving rise to a linear Stark shift. 
Rotational constants A = 10cm~^ and C = 1 cm~^ were used, with B = 1 . . . 10cm~^ varying linearly with k. 



Dipole moment component 


Coupled symmetry species 


Ma 


A'^ Ba Bb ^ Be 




A^ Bt Ba Be 




A^ Be Ba<^ Bt 



TABLE II; Couplings between the different symmetry species 
of the asymmetric rotor belonging to the four group D2 in- 
duced by the different dipole moment components, when an 
external electric field is applied along the Z-axis. 



between the difFercnt symmetry species are given in Ta- 
ble |TT] for an external electric field along the Z-axis. Fig- 
ure [5] shows the energy level structure of an asymmetric 
top molecule as a function of the asymmetry parame- 
ter K for fixed J quantum numbers. The symmetry of 
the rotational energy levels, represented by the symmetry 
species A, B^, Bi,, of the four group D2, is indicated. 
The couplings induced by the different dipole moment 
components fig according to Table [IT] are also indicated 



in the figure. 

As can be seen from Fig. [51 for a near-prolate asym- 
metric top, K w — 1, the energy levels with same Ka are 
near degenerate. These states are coupled by a dipole 
moment fia along the a-axis. This gives rise to linear 
Stark shifts, once the Stark interaction overcomes the 
splitting of these energy levels in the asymmetric rotor. 
Even in the most asymmetric case At = 0, the split- 
ting between the states with highest rotational energy, 
Ka = J, Kc = 0,1, decreases with increasing J, as can 
be seen by comparing the energy level curves for increas- 
ing values of the J quantum number. This diminishing 
splitting gives rise to linear Stark shifts in the limit of 
large J for a dipole moment Ha along the a-axis, even in 
this most asymmetric case. 

Similar arguments hold for a near-oblate asymmetric 
top, K sa 1. Here, the energy levels with same Kc are near 
degenerate. These energy levels are coupled by a dipole 
moment /ic along the c-axis, leading to linear Stark shifts 
as discussed for the near-prolate case. Once again, in the 
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most asymmetric case k = linear Stark shifts will occur 
in the limit of large J, however for different rotational 
states. The states Ka = 0, 1, Kc = J, being coupled by 
a dipole moment /_tc along the c-axis, show a decreasing 
splitting with increasing J, which can then give rise to 
linear Stark shifts. 

The situation is, however, completely different for a 
dipole moment fif, oriented along the 6-axis. As can be 
seen from Fig. [5l the energy levels which are coupled by 
such a dipole moment are always separated by another 
energy level in between. Therefore, the splitting between 
states coupled by a dipole moment along the 6-axis never 
approaches zero within one J system, independent of the 
value of the asymmetry parameter k. More severely, the 
splitting between these energy levels coupled by a dipole 
moment fib along the 6-axis even increases with increas- 
ing J. As can be seen from perturbation theory, this 
reduces the Stark shift for increasing rotational energy 
(large J). Of course, one could consider electric fields 
large enough to cause the Stark interaction to overcome 
this splitting of rotational states. Nonetheless, this gen- 
erally does not lead to linear Stark shifts for low-field- 
seeking states. Coupling to other J states becomes rele- 
vant first, forcing all states to become high- field-seeking. 

The Stark shift properties examined in the preceding 
paragraphs can also be understood from a purely classical 
point of view. For a classical rotating body, the rotation 
is stable around the axis of least and largest moment of 
inertia, i.e. the a- and c-axis. Aroimd the axis of interme- 
diate moment of inertia no stable rotation is possible. If 
the dipole moment is oriented along the a- or c-axis, as is 
the case for a symmetric top or e.g. for an a-type asym- 
metric rotor such as formaldehyde (H2CO), there exists a 
non-zero expectation value of the projection of the dipole 
moment on the space-fixed electric field axis, except for 
a classical motion with the axis of rotation perpendicular 
to the external field axis. This projection of the dipole 
moment can directly interact with the applied electric 
field, giving rise to a linear Stark shift. Conversely, if 
the dipole moment is oriented along the 6-axis, i.e. the 
axis of intermediate moment of inertia, the expectation 
value of the projection of the dipole moment on the elec- 
tric field axis will be zero. Hence, the external electric 
field first has to hinder the rotation and orient the dipole, 
which becomes more and more difficult with increasing 
rotational energy and angular momentum. Therefore, 
the magnitude of the oriented dipole decreases with in- 
creasing rotational energy, in accordance with the de- 
scription given by quantum mechanics. This oriented 
dipole moment can then interact with the electric field. 
This second-order interaction gives rise to a quadratic 
Stark shift for 6- type rotors. The same argument holds 
for linear heteronuclear molecules with a permanent elec- 
tric dipole moment, where the rotation averages out the 
projection of the dipole moment. Hence, linear molecules 
exhibit a quadratic Stark effect unless an electronic angu- 
lar momentum parallel to the electric dipole moment /i is 
present, as e.g. in H states of NH, OH or CO*. However, 
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TABLE III: Nuclear spin degeneracy factors gi for the three 
water isotopologs H2O, D2O and HDO [s^]. e and o refer to 
the parity (even or odd) of the quantum numbers Ka and Kc 
respectively. 

note that in the case of accidental degeneracies also a b- 
type asymmetric rotor may show linear Stark shifts for 
certain states [1^ . The mechanical analogue of these ac- 
cidental degeneracies, if existing, is not evident. Overall 
one can summarize that even molecules such as the three 
water isotopologs H2O, D2O and HDO which seem very 
similar at a first glance, can show surprisingly different 
behavior when exposed to an external electric field. 

IV. CALCULATION OF THE GUIDED FLUX 

With the Stark shifts of the different isotopologs at 
hand we are now ready to calculate their guided fluxes. 
The total guided flux of a molecular species can be cal- 
culated by summing over the contributions to the flux of 
all individual internal states, weighted with their thermal 
occupation in the source. The thermal occupation of a 
rotational state | J, r, M) in the source is given by 

PJtM = -^9m9i exp {-EjrM /ksT), (10) 
with the partition function 

Z = ^ 9Mgi exp {-EjrAi/kBT). 

JtM 

EjtM is the rotational energy, T the source temperature, 
gM the M degeneracy factor of the state, and gi the 
nuclear spin degeneracy factor. Nuclear spin degeneracy 
factors for the different isotopologs are listed in Table Hill 

The guided flux for a given guiding electric field can 
now be readily expressed as 

$ = A^O ^ PJrAlfiAWjrM), (H) 
JtAI 

with the guidable fraction of a rotational state \JtM) 
given by f{AWj^nj) oc (AVFj^j^^)^ as shown in Sec- 
tion [ni and Nq the number of molecules injected into 
the guide per second. Figure [6] compares Stark shifts 
and calculated source populations of individual rota- 
tional states for H2O and D2O. From this figure it can 
already be anticipated that in the case of H2O the state 
I J = 1, r = 1, M = 1) will contribute strongly to the 
guided flux. The reason for the large thermal population 
of this state as compared to the | J = 1, r = 1, M = 1) 
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FIG. 6: (Color online). Calculated thermal populations in 
the source PjtM at a source temperature of 293 K, and Stark 
shifts AWj^M ^-t Bill electric field of 130kV/cm for the ro- 
tational states \ J,T,M) of H2O and D2O. States with large 
contributions to the guided flux are labeled with quantum 
numbers. 



state in D2O can be found in the different nuclear 
spin statistics of the two isotopologs. For such consid- 
erations based on symmetry properties it is advanta- 
geous to transform from | J, r, M) to the quantum num- 
bers I J, Ka, Kc, M). For the state | J = 1, r = 1, M = 1) 
this results in the corresponding quantum numbers 
\J =l,Ka^l,Kc = Q,M =1)- According to Table Hm 
this state is favored in H2O by 3 : 1, while it is unfavored 
in D2O by 3 : 6. 

The expected guided fluxes of the different isotopologs 
calculated this way are depicted in Fig. [7] as a func- 
tion of the electric field. The guided fluxes of the dif- 
ferent isotopologs differ remarkably. Furthermore, not 
only the flux but also the electric field dependence is 
clearly different. While H2O and D2O show a quar- 
tic dependence of the guided flux, caused by quadratic 
Stark shifts, the electric field dependence of HDO shows 
a mainly quadratic behavior, indicating the presence of 
linear Stark shifts. 

The calculations of the guided fiuxes allow us to de- 
duce the populations of individual rotational states in 
the guided beam. In a recent experiment these cal- 
culated state populations were experimentally verified 
by coUinear ultraviolet spectroscopy in a cold guided 
beam of formaldehyde H2CO [131 . For the water iso- 
topologs, populations of individual rotational states with 
contributions larger than 5% of the total flux of the 
considered isotopolog are listed in Table IIVI Remark- 
ably, the single state | J = 1, r = 1, M = 1) of H2O con- 
tributes «80% to the guided flux using a room tem- 
perature source. Similarly, in D2O the most populated 
state I J = 1, r = 1, il/ = 1) contributes w21 % to the to- 
tal flux. The reason for these large populations as com- 
pared to a molecule such as formaldehyde (H2CO) [l^l 
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FIG. 7: (Color online), a) Calculated guided fluxes of the 
different water isotopologs as a function of the trapping elec- 
tric field. The symbols are the results of the calculation de- 
scribed in Section ITVl whereas the solid curves are quartic and 
quadratic fits to the calculated flux. Guided fluxes are nor- 
malized to the flux of H2O at a trapping field of lOOkV/cm. 
The guided flux of HDO shows a quadratic dependence on the 
trapping electric field, whereas for H2O and D2O a quartic de- 
pendence is found, b) Zoom-in to compare the calculations 
of the guided flux for H2O and D2O to a quartic curve. 



with similar size rotational constants but mainly linear 
Stark shifts can be found in the quadratic Stark shift be- 
havior. As can be seen from Fig. 21 the size of the Stark 
shift decreases with rotational energy for a molecule with 
a quadratic Stark shift, leaving only few low-energy states 
with largest Stark shifts for guiding psj . 



V. EXPERIMENTAL SETUP 

To compare with theory, experiments with the three 
water isotopologs H2O, HDO and D2O were performed. 
As shown in Fig. El a setup similar to the one described 
previously [2^, l27f is used for the experiment. Water 
molecules are injected into the quadrupole guide through 
a ceramic tube with diameter 1.5 mm and a length of 
9.5 mm. The pressure in the reservoir is kept at a fixed 
value of 0.10 mbar via a stabilized fiow valve, resulting 
in a gas flow of 1 x 10~^ mbar-l/s. Since water has a 
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J T M 


Erot PS AW PG 

(cm-1) (%) (cm-i) (%) 


1 1 1 

2 2 

3 1 3 


42.3 2.9 0.13 79.3 
95.7 0.8 0.08 7.4 
213.6 1.3 0.05 5.1 


D2O 


J T M 


Erot PS AW" PG 

(cm-1) (%) (cm-i) (%) 


1 1 1 

3-2 1 

2 2 
3-2 
3-2 2 


22.7 0.54 0.20 21.3 
74.5 0.83 0.16 21.2 
49.3 0.94 0.13 17.1 
74.5 0.42 0.18 13.2 
74.5 0.83 0.10 8.0 


HDO 


J T M 


Erot PS AVK" PG 
(cm-i) (%) (cm-i) (%) 


3 3 3 

2 2 2 

4 4 4 

3 3 2 

4 2 4 


234.9 0.45 0.83 15.0 

109.8 0.83 0.59 14.3 
407.0 0.19 0.89 7.5 

234.9 0.45 0.55 6.6 
297.8 0.33 0.58 5.5 



TABLE IV: Properties of selected rotational states (popula- 
tion in the guide>5 %) of the three water isotopologs H2O, 
D2O and HDO. Erot' zero field rotational energy, pg: Source 
population for a temperature of 293 K, AW: Stark shift at 
an electric field of lOOkV/cm, pG'- population in the electric 
guide at an apphed electric field of lOOkV/cm. 
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FIG. 8: (Color online). Experimental setup, a) Water 
molecules from the thermal reservoir enter the electric guide. 
Slow molecules are trapped in the quadrupole field and guided 
to an ultra-high- vacuum chamber, where they are detected 
by the quadrupole mass spectrometer, b) Electric field distri- 
bution in the quadrupole guide for ±5kV electrode voltage, 
resulting in a trapping electric field exceeding 93kV/cm. 



sufficiently high vapor pressure of ~ 25 mbar at room 
temperature [37| , no heating of its container is necessary. 
The constituents of the gas injected through the tube 
can be monitored by a residual gas analyzer placed in 
the source vacuum chamber. The guide is composed of 
four stainless steel electrodes of 2 mm diameter separated 
by a distance of 1 mm. For an electrode voltage of ±5 kV, 
a trapping electric field of around 93 kV/cm is generated. 
Transversely slow molecules in a Ifs state are trapped by 
the enclosing high electric fields. A longitudinal cut-off 
velocity is obtained by bending the guide. Slow molecules 
are guided around two bends with a radius of curvature 
of 5 cm and through two differential pumping stages to an 
ultrahigh-vacuum chamber, where they are detected by 
a quadrupole mass spectrometer (QMS). The molecules 
arc ionized by electron impact in a cross-beam ion source 
and mass selected in the analyzer. In the final stage of 
the QMS single ion counting is performed. 

Measurements of D2O were performed with isotopi- 
cally pure D2O. For measurements of HDO, mixtures of 
H2O and D2O with a ratio of 1:1 and 4:1 were used. 
These resuh in a HDO fraction of «48% and «27%, re- 
spectively, coming from the source as measured with the 
residual gas analyzer. When using these mixtures in the 
source, all three isotopologs are guided simultaneously. 
However, due to its much larger Stark shift for some ro- 
tational states (see Section iHll and Fig. |4]), HDO is pref- 
erentially guided. This measurement with the mixture 
of different isotopologs allows to extract contributions of 
HDO and D2O. The contribution of H2O is also visible, 
although it is shadowed by fragments of D2O and HDO 
in the QMS at mass 18, 17 and 16 amu where H2O is de- 
tected. Measurements of H2O were performed with pure 
H2O to avoid these unwanted contributions of D2O and 
HDO. 



VI. RESULTS AND DISCUSSION 

The experiments are performed as a series of time of 
flight measurements. Here, the high voltage (HV) is 
switched on and off repeatedly in a fixed timing sequence. 
A typical time of flight trace is shown in Fig. [9l To sub- 
tract background contributions, the guided signal for the 
different isotopologs is determined from the difference in 
the steady state QMS signal with HV applied to the guide 
and HV switched off. 

Before comparing the signals of guided molecules to 
the predictions of the theory presented in Section HVl the 
detection process is considered in some more detail. The 
QMS used for detection of the guided molecules operates 
as a residual gas analyzer, which measures the densities 
of individual constituents in the recipient. The signals of 
guided molecules observed in previous [l^, HI] as well as 
in the present experiments are, however, fully compatible 
with the predicted electrode voltage dependence of the 
guided flux. 

The solution to this surprising fact lies in the pros- 
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FIG. 9: (Color online), a) Time of flight trace measured with 
HDO for an electrode voltage of ±5 kV. High voltage (HV) is 
switched on for an interval of 250 ms, with a 50 % duty cycle. 
The QMS signal starts rising with a time delay after switching 
on HV, until a steady state is reached. This time delay cor- 
responds to the fastest molecules' travel time from the nozzle 
to the QMS. The steady-state value is used for voltage de- 
pendence measurements. The slow decay after switching off 
the high voltage is due to the fact that the last section of the 
guide closest to the QMS is not switched to avoid electronic 
pick up on the QMS. Hence, molecules in this part can arrive 
at the QMS even after the HV has been switched off in the 
first section, b) Zoom into the rising slope of the signal, from 
which a velocity distribution can be derived. These plots show 
raw data, i.e. no background subtraction has been applied, to 
illustrate the excellent signal to noise ratio. 



sure dependence. The model of velocity filtering pre- 
sented in Section [ill and IIVI assumes that the molecules 
from a thermal ensemble are transferred into the guide 
while preserving their velocity distributions. Measure- 
ments performed with deuterated ammonia ND3 for very 
small inlet pressures indeed show the characteristic volt- 
age dependence expected for a density measurement [ssj . 
For increasing inlet pressure, however, the voltage depen- 
dence changes, and the signal more and more resembles 
a flux measurement. Here, collisions of molecules in the 
nozzle and in the "high pressure" region directly behind 
the nozzle become more likely. Slow molecules are more 
likely removed from the ensemble, effectively leading to 
a "boosting" of the beam. This boosting is also visible in 
the velocity distributions of the guided beam presented 



in Section Em 

From the pressure dependence observed in these mea- 
surements with ND3 the following conclusion can be 
drawn: by including this boosting effect in the model, 
the voltage dependencies over the entire pressure range 
studied agree with a density measurement. For the con- 
ditions chosen in previous experiments as well as in the 
experiments described in this paper the electrode voltage 
dependence of the QMS signal resembles a flux measure- 
ment of a gas without boosting. Therefore, the outcome 
of the experiment can be directly compared to the theory 
for an ideal effusive source presented in Section HIl and [TVl 
To avoid any ambiguity, we refer to the signal of guided 
molecules as " molecule signal" or " detector signal" . 

To extract the relative molecule signals of the differ- 
ent isotopologs from the individual measurements, some 
corrections are necessary. First of all, the fragmenta- 
tion of the guided molecules in the QMS has to be taken 
into account, since the measurement is performed at only 
one mass (18amu for H2O, 20amu for D2O, and 19amu 
for HDO). The correction is determined from the ion 
count rates of the different fragmentation products for 
pure H2O (18, 17, 16amu), pure D2O (20, 18, 16amu), 
and for the mixture with high (48%) HDO content (19, 
18, 17, 16amu), where HDO dominates the molecule sig- 
nal. Secondly, the ion count rate for each isotopolog is 
corrected for the relative contributions of the isotopologs 
injected into the guide as monitored by the residual gas 
analyzer in the source chamber. A velocity distribution of 
the guided molecules is determined from the rising slope 
of the time of flight trace. 



A. Electrode Voltage Dependency 

As shown in Fig. [101 the electrode voltage dependence 
of the molecule signals for H2O and D2O are well de- 
scribed by the calculations over more than two orders 
of magnitude. The calculations are scaled by only one 
global scaling factor to account for detection efficiencies 
and the amount of gas injected through the nozzle into 
the guide. Note that no relative scaling factor between 
the calculations for H2O, D2O and HDO has been em- 
ployed. This excellent agreement between experiment 
and theory verifies that the filtering process is well de- 
scribed by the model presented in Sections [Til HIH and 

nvi 

As can be seen from Fig. [TUl HDO shows the largest 
molecule signal of all three water isotopologs, as pre- 
dicted from calculations. We experimentally determine 
a ratio of detector signals between HDO and D2O at a 
guiding field of 130 (100) kV/cm of 16.4 (18.9), where the 
calculations predict 14.7 (20.5). This is well within the 
overall uncertainty given e.g. by the determination of the 
relative contributions of isotopologs injected through the 
nozzle as derived from the residual gas analyzer signal. 
The measurements done with different HDO amounts 
give count rates agreeing to within 5-10%, which allows 
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FIG. 10: (Color online). Guided signal of H2O, D2O and 
HDO as a function of the applied electrode voltage. Shown 
as solid curves are the theoretically predicted signal depen- 
dencies, adjusted to fit D2O data at ±6kV electrode voltage. 
Note that only one global scaling factor is used for all theory 
curves. The different slopes of the curves in the double log- 
arithmic plot directly indicate the different Stark shift prop- 
erties. H2O and D2O show a quartic dependence of guided 
signal on the applied electric field, indicating quadratic Stark 
shifts, whereas for HDO a mainly quadratic dependence indi- 
cating linear Stark shifts is found. 



an estimate for the uncertainty of the residual gas ana- 
lyzer corrections applied. Furthermore, the guided flux 
of HDO largely stems from rotational states at higher ro- 
tational energies. These are therefore more strongly af- 
fected by centrifugal distortion corrections not included 
in the rigid rotor approximation being used, leading to 
energy level shifts and hence changes of Stark shifts which 
can affect the accuracy of the calculations. 

The measured signals of guided H2O and D2O show 
a quartic dependence on the applied electrode voltage, 
as was shown for D2O already in a previous experiment 
[2^. The electrode voltage dependence for the molecule 
signal of HDO is best described by a quadratic behavior, 
confirming the main contributions from states with lin- 
ear Stark shifts predicted by calculations. This different 
dependence of the guided signal on the applied electric 
field is directly evident from different slopes in the dou- 
ble logarithmic plot shown in Fig. [TOl Also the good 
agreement between calculations and the experiment over 
a wide range of applied electrode voltages resulting in 
changes of several orders of magnitude for the molecule 
signal is remarkable. 

The detector signals of the water isotopologs can also 
be roughly compared to other molecules used so far. For 
deuterated ammonia ND3, count rates of the order of 
3 X 10^ cts/s were observed in the same setup at an elec- 
trode voltage of ±5kV, using a reservoir at room tem- 
perature. This is to be compared to a guided signal of 
w 4 X 10* cts/s for HDO. The reason for this smaller flux 
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FIG. 11: (Color online). Normalized velocity distributions for 
the water isotopologs H2O, D2O and HDO at an electrode 
voltage of ±5kV. The velocity distributions are derived from 
the rising slope of time of flight measurements. The solid 
curves are a guide to the eye. Vertical arrows indicate the cut- 
off velocities obtained by linear extrapolation (dashed lines) 
of the high velocity side. 



can be found in the fact that for HDO only few states 
exhibit linear Stark shifts and hence contribute to the 
guided flux, whereas ND3 exhibits mainly linear Stark 
shifts. Furthermore the Stark shifts of HDO arc smaller 
due to the smaller dipole moment component along the 
a-axis as compared to the dipole moment of ND3. Nev- 
ertheless, the fact that molecules with so different Stark 
shift properties can be efficiently filtered out of a ther- 
mal gas illustrates the wide applicability of the velocity 
filtering technique. For chemically stable polar molecules 
which can be supplied to the nozzle via the teflon tube, 
it is sufficient to bring them into the gas phase in a reser- 
voir at a pressure in the 0.1-1 mbar range. At smaller 
pressures guiding is still possible, however at the expense 
of reduced signals. Chemically very reactive molecules 
such as molecular radicals (typical examples are OH and 
NH) arc in principle also guidable. Beams of molecu- 
lar radicals can be loaded into a helium buffer gas cell, 
where they thcrmalizc by collisions with the cryogenic he- 
lium gas [3^. The translationally and internally cooled 
down molecules can then be extracted by the electrostatic 
quadrupole guide (29j and made available for further ex- 
periments. 



B. Velocity Distributions 

The different Stark shift properties of the various iso- 
topologs are also evident from their velocity distribu- 
tions. The cut-off velocity, i.e. the velocity of the fastest 
molecules which can still be guided, is given by the 
molecules' Stark shift. Hence, it depends on the rota- 
tional state of the molecule. Since for molecules from 
a thermal source many states contribute to the guided 
flux, the cut-off velocity is given by the state with the 
largest Stark shift. In the experiment the cut-off velocity 
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is determined from a linear extrapolation of the high- 
velocity side of the velocity distribution towards zero. 
From the measurements performed at an electrode volt- 
age of ±5kV it can be seen in Fig. llll that the maximum 
of the distribution as well as the cut-off velocity shifts 
towards higher velocities from H2O and D2O to HDO. 
This is caused by the much larger Stark shifts of HDO 
as compared to H2O and D2O (see Fig. S]). The wide ve- 
locity distribution for HDO is caused by the large Stark 
shifts found for this isotopolog and by the used bend ra- 
dius of 5 cm. We experimentally determine cut-off veloc- 
ities of 60m/s for H2O, 69m/s for D2O, and 130m/s for 
HDO at an electrode voltage of ±5 kV. Calculated Stark 
shifts of states contributing to the guided flux (similar 
to Table Hvl but for a guiding field of 93kV/cm reached 
at ±5 kV electrode voltage) result in cut-off velocities of 
57m/s for H2O, 67m/s for D2O and 125 m/s for HDO, 
which is in good agreement with measurement. For HDO 
states with large Stark shifts of up to 0.80 cm^^ at ±5 kV 
electrode voltage are predicted to contribute with « 20% 
to the guided flux. These states might be responsible for 
the small but non-zero signal in the HDO velocity distri- 
bution between 130m/s and 150m/s, being supported by 
the fact that a cut-off velocity of 150 m/s corresponds to a 
Stark shift of 0.80 cni~^. Note that this signal exceeding 
the linear extrapolation to the falling slope of the veloc- 
ity distribution is present in the velocity distributions of 
HDO measured at ±3kV, ±5kV and ±7kV. 

On the low velocity side of the velocity distribution 
a linear extrapolation does not cut the horizontal axis 
at zero velocity. We attribute this to collisions of the 
molecules in the vicinity of the effusive source, which 
removes the slowest molecules from the thermal distribu- 
tion, effectively leading to a "boosting" of the molecules 
entering the guide. This is supported by measurements 
with other gases such as ND3 , where an influence of the 
reservoir pressure on the number of molecules with low 
velocities was observed. 

Measurements of velocity distributions for different 
electrode voltages (±3kV, ±5kV, ±7kV) were per- 
formed for D2O and HDO, and are shown in Fig. [T21 
In the data a shift of the maximum of the distribution 
and of the cut-off velocity, i.e. the maximum velocity of 
molecules which can still be guided, can be seen. This 
is to be expected, since larger voltages applied to the 
guide electrodes result in a deeper trapping potential, 
hence faster molecules are guided. Regarding the depen- 
dence of the cut-off velocities on the applied electrode 
voltage in Fig. [T31 a linear dependence in the case of 
D2O and a square-root dependence in the case of HDO 
is observed. As was shown in Section [Hi the cut-off ve- 
locity depends on the square-root of the Stark shift at 
the applied guiding fleld. Hence, a linear dependence as 
for D2O is found for molecules with quadratic Stark shift, 
while the square-root dependence is found for molecules 
as HDO where states with linear Stark shifts dominate 
the guided flux. 
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FIG. 12: (Color online). Normalized velocity distributions 
for D2O and HDO as a function of electrode voltage (±3kV, 
±5kV, ±7kV). The solid curves are a guide to the eye. The 
vertical arrows indicate the cut-off velocity, which is deter- 
mined from a linear approximation (dashed lines) to the high 
velocity side. 




Electrode voltage (kV) 

FIG. 13: (Color online). Voltage dependence of cut-off ve- 
locities for D2O and HDO. For D2O the line is a linear fit, 
for HDO the solid curve is a fit of a square root dependence. 
These dependencies are expected for quadratic (D2O) and 
linear (HDO) Stark shifts respectively. 
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TABLE V; The three possible connections between the prin- 
cipal axes of the molecular moment of inertia tensor (a, b, c) 
and the right-handed body-fixed coordinate system {x, y, z). 
The ^ denotes the choice of right-handed coordinate systems 
(adapted from [sF 
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TABLE VI: Coefficients used in the matrix elements of the 
reduced Hamiltonian (Eq. O for the different representations 
listed in Table |V](adapted from [H). 



VII. SUMMARY AND OUTLOOK 

To summarize, we have produced continuous cold 
guided beams of the water isotopologs H2O, D2O and 
HDO by electrostatic velocity filtering. We discuss in 
detail the influence of molecular parameters such as ro- 
tational constants and orientations of electric dipole mo- 
ments on the behavior of the molecule in an external 
electric field. Based on this discussion of molecular Stark 
shifts, we develop the theory for the electrostatic veloc- 
ity filtering process. Wc find excellent agreement be- 
tween the experimentally observed signals and our cal- 
culations of guided fluxes, confirming that the filtering 
process is well described by the presented model. The 
molecule signals of the different isotopologs show quar- 
tic and quadratic electrode voltage dependencies respec- 
tively caused by quadratic Stark shifts for H2O and D2O 



and by linear Stark shifts of the states contributing most 
to the guided flux in HDO. These different Stark shift 
properties can also be seen from the dependence of the 
cut-off velocity on the applied guiding electric field as 
measured for D2O and HDO. Furthermore, calculations 
of guided fluxes allow to deduce populations of individ- 
ual rotational states in the guided beam. These can be as 
high as 80% in the case of H2O for a room-temperature 
source. Overall, this shows that velocity filtering by an 
electric guide is a technique applicable to a wide variety 
of molecular species. 

In the presented experiment, the internal state dis- 
tribution was only inferred from calculations based on 
the good agreement between theory and experiment. 
However, internal state diagnostics of cold guided water 
beams or any other guided species should be possible by 
transferring the depletion spectroscopy technique used 
for formaldehyde in the near-ultraviolet [27[ to either 
the infrared using vibrational transitions or to the mi- 
crowave domain using purely rotational transitions. Such 
additional state-dependent detection will be beneficial for 
studies of, e.g., collision processes with or between cold 
molecules. Furthermore, with cold water beams at hand 
it should be possible to investigate in the laboratory pro- 
cesses such as, e.g., ice formation on dust grains under 
conditions present in the interstellar medium. 
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APPENDIX: TABLES 

In this appendix we provide tables IV^IVHI 



J 





J+l 


J' 

J 


J- 1 


(^1 \J'} 

{JK\ iS>F, \J'K) 

{JK\^Fy\J\K±l) 

= ^i{JK\^F^\J',K±l) 
{JM\^za\J'M) 

(JM|<I>yg|J',Af±l> 
= ±i{JM\ ^Xg \J',M±1) 


[4(J + 1) [(2J + l)(2J + 3)]'/']~' 
2[(J + 1)^-A-^]'/' 
=F[(J±A'+1)(J±A + 2)]^/" 

=F[(J±A/+1)(J±A/ + 2)]'''^ 


[4J(J + 1)]-^ 
2A' 

[(JtA') (J±A-f 1)]^/^ 
2M 

[(J=p Af) (J± Af + l)]'/^ 


-2 [J^-A^]^/' 
T[(JtA') (JtA'-1)]^/^ 
-2 [J2 - A/2] 

t[(Jtm) (Jtm-2)]^/2 



TABLE VII; Direction cosine matrix elements {J,K,M\ $Fg \ J'K'M'} for the symmetric rotor (adapted from [s^). 
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